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The effect of disorder on the Landau ieveis of massless Dirac fermions is examined for the cases 
with and without the fermion doubling. To tune the doubling a tight-binding model having a 
complex transfer integral is adopted to shift the energies of two Dirac cones, which is theoretically 
£\J ■ proposed earlier and realizable in cold atoms in an optical lattice. In the absence of the fermion 

doubling, the n — Landau level is shown to exhibit an anomalous sharpness even if the disorder is 
uncorrelated in space (i.e., large K-K' scattering). This anomaly occurs when the disorder respects 
the chiral symmetry of the Dirac cone. 
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Kicked off by the experimental observation of the graphene quantum Hall effect 1 ' 2 , fascination with massless Dirac 
fermions is mounting, where they appear not only in graphene but more generically in various systems such as organic 
■ metals^r— , cold atom systems in optical lattices^ and molecular graphene^. Among these systems, the number of 
| massless Dirac fermions is always even for solid state materials, which is called the "fermion doubling"-'-. On the 
other hand, manipulation of Dirac cones into single cones has been theoretically considered^, which may be realizable 
in optical lattices where the Hall conductivity is detectable experimentally^. 

Inspired by these, we explore here the effect of disorder for the Landau levels of massless Dirac fermions, in particular 
in the absence of the fermion doubling. Specifically, the effect on the n = Landau level, which is essential to the 
anomalous Hall effect of massless Dirac fermions, is examined from the viewpoint of the fermion doubling and the 
symmetry of the system. In our previous worki^, we have shown that the n = Landau level of the honeycomb lattice 
(graphene) becomes anomalously sharp even in the presence of disorder, if the disorder respects the chiral symmetry 
and is spatially correlated over a distance exceeding a few lattice constants. Conversely, we have also pointed out that 
when the disorder is spatially uncorrelated, the n = Landau level is broadened just like the other Landau levels, 
even if the chiral symmetry is respected by the disorder. 

To explore the single-to-double Dirac cone crossover, here we consider a two-dimensional lattice model having two 
Dirac cones which are shifted in energy with each other as in the model proposed by Watanabe et alJ^. We examine 
numerically the effect of disorder with this lattice model, and have found that, even if the disorder is uncorrelated in 
space, the n = Landau levels start to become anomalously sharp as the two Dirac cones are energetically shifted. 
Notably, the anomalous sharpness in the absence of the fermion doubling occurs when the disorder respects the chiral 
symmetry for each cone. In fact, a potential disorder, which breaks the chiral symmetry, washes out the sharpness. 
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II. MODEL AND NUMERICAL RESULTS 



We adopt here a two-dimensional square lattice with the following nearest-neighbor (NN) t and the next nearest- 
neighbor (NNN) t' transfer integrals, 

X 

d ' H = — £c r+e ^c r + (— 1) A+ y tc rJre ^c r 

r 

+it'{ct +ex+ey c r ) + R.c, t,t'eR (1) 

where r — (n x ,n y ) denotes the lattice points and e x = (l,0)(e y = (0, 1)) the unit vector in x(y) direction with all 
lengths measured in units of the lattice constant. To realize shifted Dirac cones some transfer energies have to be 
complex, and here the NNN transfer is pure imaginary. Although complex transfer integrals may seem unrealistic, 
they can be realized in cold atoms in optical lattice d 11 ' 13 . In the absence of a magnetic field, the Hamiltonian in the 
momentum space is expressed as 
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FIG. 1: Left: The present lattice model. Right: Energy dispersion in the fci-fe plane for t' = OAt. 




FIG. 2: Density of states as function of the Fermi energy E for various values of i'(oc the energy shift of two Dirac cones) 
averaged over 5000 samples with a system size 30 by 30. The magnetic flux piercing each square plaquette is 0.01(/i/e) and 
the strength of disorder a/t = 0.1. Inset: Density of states for the same parameters when we replace the bond disorder with a 
potential disorder with a variance O.lt. 

where A(fc) = — 1(— 1 + e tkl + e lkl+ 2 + e lk2 ) with k\ = k ■ e.\ and fc 2 = k ■ e 2 with the primitive vectors taken to 
be ei = e x — e y and e 2 = + e y . In this model, we have two Dirac cones at (fci, fc 2 ) = i(7r/2, —tt/2) with energies 
±2t' , so that the two Dirac cones are shifted in energy from each other when the strength t' is non-zero. A magnetic 
field is introduced by taking the Peierls substitution t{t') — > t(t')e~ 27!ie( > r \ where the summation of the phases 9{r) 
along a closed loop is equal to the magnetic flux enclosed by it in units of h/e . We introduce a random component 
6t(r) for the NN transfer integral i(r) = t + St(r) that has a gaussian distribution with variance a and is uncorrelated 
in space (St(r)St(r')) = <r 2 6(r — r'). Still, this disorder preserves the chiral symmetry for each Dirac cone. The bond 
disorder considered in our previous work 12 is equivalent to the present disorder in the limit of zero correlation length. 

The numerically evaluated density of states in a magnetic field is shown in FigJ^l where both n = and n = ±1 
Landau levels are split as the two Dirac cones are energetically shifted with the increase of t' . Remarkably, the n = 
(zero-mode) Landau level for each cone becomes anomalously sharper with t' even in the presence of disorder, while 
the n = ±1 Landau levels remain broadened. For comparison, the density of states calculated for a potential disorder, 
which breaks the chiral symmetry, shows no such anomaly (FigJ^l inset). 

III. SUMMARY AND DISCUSSIONS 

We have numerically investigated the effect of disorder on the Landau levels of the massless Dirac fermions with 
and without the fermion doubling. We have clearly shown that, if the chiral symmetry for Dirac cones is respected, 
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the zero energy (n — 0) Landau level becomes anomalously sharp in the absence of the fermion doubling even when 
the disorder is uncorrelated in space. Since the shift in energy of the Dirac cones suppresses the mixing between the 
Dirac points due to the disorder scattering, the present result implies that the broadening reported for uncorrelated 
disorder in the honeycomb lattice (graphene)i2 is due to the mixing between the two Dirac fermions with opposite 
chirality. 
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